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Abstract. We use the exact determinantal representation derived by Kitanine, Maillet, and Terras for ma- 
trix elements of local spin operators between Bethe wave functions of the one-dimensional s = i Heisenberg 
model to calculate and numerically evaluate transition rates pertaining to dynamic spin structure factors. 
For real solutions z\, . . . , z T of the Bethe ansatz equations, the size of the determinants is of order rxr. We 
present applications to the zero-temperature spin fluctuations parallel and perpendicular to an external 
magnetic field. 



1 Introduction 

The one-dimensional (ID) s = | Heisenberg model was 
solved, in principle, by Bethe some 70 years ago using 
an ad-hoc trial wave function, now famously known as 
the Bethe ansatz Decades went by before the power 
and scope of this method of exact analysis became widely 
known for applications to spectrum and thermodynamics 
of a select class of completely integrable model systems 
§■ 

Until recently, the most vexing exception to immense 
progress in the further development of the Bethe ansatz 
has been the absence of a practical method to use the 
exactly known and readily available Bethe wave functions 
for the explicit calculation of transition matrix elements. 
The knowledge of such transition rates is of paramount 
importance for an understanding of dynamic correlation 
functions in relation to the underlying quasiparticles and 
for the interpretation of experimental probes of quantum 
fluctuations in quasi-lD magnetic compounds. 

It was most remarkable, therefore, when Kitanine, Mail- 
let, and Terras || succeeded in reducing matrix elements 
between Bethe wave functions for local spin operators to 
determinantal expressions. Here we use these expressions 
and the norms previously determined by Korepin H to 
calculate dynamic spin structure factors of the Heisenberg 
antiferromagnet with periodic boundary conditions and a 
magnetic held: 
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2 Transition rates expressed as determinants 

Consider Bethe wave functions with z-component Sj. — 
N/2 — r of the total spin. They are specified by sets of ra- 



pidities Z\, . . . ,z r , which are solutions of the Bethe ansatz 
equations 

N<t>{z i ) = 2iTli + ^2<l>[{zi-z j )/2], i = l,...,r, (2) 

where 4>(z) = 2arctanz. The Bethe quantum numbers Ii 
provide a natural classification of the spectrum. The chal- 
lenge is to calculate transition rates for physically moti- 
vated operators from the information encoded in the Zi. 
Using the Bethe wave function directly is feasible but com- 
putationally inefficient and therefore limited to relatively 
small systems {r< 12) ||. The combined advances re- 
ported in Refs. are the basis of a much more powerful 
approach, which will be used in the following to calculate 
transition rates 
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from the ground state of (0) for the operators 

Sq = E/ eiqn ^n> A* = 2>+j ~ ■ 



(3) 



(4) 



They probe the parallel (// = z) and the perpendicular 
(/i = +, — ) spin fluctuations at zero temperature. 

The transition rates for the parallel spin fluctuations as 
inferred from Eq. (5.12) of Ref. || for the matrix element 
(^oIS^I^a) can be brought into the form 
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where 
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Rescaling H, P consolidates powers of two. The calculation 
of the transition rate is thus reduced to the evaluation of 
three r x r determinants: 
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£ r ({^})=n K (^)' ^r({^})=n^-^-) 5 (17) 
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P ab = i2n(z° a ) J] G{z 3 - z h ), 



(19) 



The rapidities {z®} and {z.;} belong to the ground state 
l^o) and to the excited state \ipx), respectively. The norm 
of \ip\) from Ref. M, transcribed to our notation, reads: 



G(z) = - + i. 



(20) 
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Along similar lines of manipulation, the transition rates 
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and can be compactified into 

2 r2 dctK({zJ) 
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where 
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for the perpendicular spin fluctuations as inferred from 
the Eqs. (5.3) and (5.7) of Ref. || for the matrix element 
(•00 1 $n IV'a) are brought into the form: 
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For real 2, the absolute square of (pi) becomes: 



(14) 



H+ _ +1 = in(z a ), a= l,...r + l, b= 1,... ,r, (23) 
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(24) 
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3 Applications 

To demonstrate the efficacy of expressions (|l^) and ( p2| ) 
for the transition rates M^(q), /i = z, +, — , we present 
three applications to the dynamic structure factors 

S0i(q,u)=2irY,M£(q)6(u-ux). (25) 

A 

We begin with the parallel spin fluctuations at magnetiza- 
tion M z /N = i (half the saturation value). The spectral 
weight of S zz (q,uj) was shown to be dominated by a set 
of collective excitations consisting of two unbound quasi- 
particles named ip and ip* ||. The spectral range of the 
ipip* continuum, shown in Fig. |l|(a), has soft modes at 
q = 0, 7r/2 and is partly folded back onto itself along a 
stretch of its upper boundary. 
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Fig. 1. (a) Energy versus wave number of the tptp* excitations 
for TV — > oo. (b) Scaled transition rates between the ground 
state and the ipip* states at q — n/2 for N = 12, 16, ... ,32 
(open circles) and N = 512 (closed circles). The dashed line is 
a fit a + buj 71 ' 2 of the N = 512 data at lo/J < 0.5. All results 
are for M z = N/4. 

In Fig. Q(b) we have plotted the scaled transition rates 
Mff*(n/2,w) = NM^(n/2) versus u pertaining to the 
ipip* excitations. The open circles are results previously 
obtained for N < 32 by calculating matrix elements di- 
rectly from Bethe wave functions H. The full circles are 
data for N = 512 as obtained from the determinantal ex- 
pression (|l6|). 

The function Mff (n/2,oj) varies smoothly across the 
entire ipip* continuum including the fold. The observed in- 
frared singularity, Mff (tt/2^ u) ~ uj r '- 2 , rj-2 = -0.468 . . . 
confirms exact predictions 0,|J. At the other end of the 
continuum, Mff (n/2, ui) tends to go to zero or a value 
close to zero. 

The density D^" (n/2, w) = 2n / [N (uj i+1 - Ui)} of ipijj* 
states is shown in Fig. ||(a). It is flat and featureless except 
near the upper band edge, where the fold in the contin- 
uum produces a square-root divergence [[|. The spectral- 
weight distribution Sf^ (n/2, cu) resulting from the prod- 
uct of Mff and is then a double-peak structure as 
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Fig. 2. (a) Density of ipip* states and (b) spectral weight distri- 
bution of the ipip* states in Szz(q,u) at q = n/2 and M z = N/4 
from data for N — 512. 

shown in Fig. ||(b) with the divergences at the lower and 
upper band edges caused by the transition rates and the 
density of states, respectively. At frequencies where the 
continuum is folded back, the spectral weight of two lines 
of excitations must be added up to produce the correct 
lincshape, causing a (barely visible) singularity within the 
band. 

Now consider the perpendicular spin fluctuations as 

described by S \.(q,ui), again at M z /N = j. At the zone 

boundary (q — n), some 98% of the spectral weight is 
carried by a continuum of collective excitations consist- 
ing of two ip quasiparticles || . The shape of the ipip con- 
tinuum is shown in Fig. |^(a) The scaled transition rates 
M^(n,uj) = NM^(n) pertaining to the tpip excitations 
are plotted versus w in Fig. ||(b) for N = 1536 by n umer- 
ical evaluation of the determinantal expression (G3) . Also 
shown are data points for N — 12, 16, ... , 28 previously 
obtained by a different method ||. 

Again we find a smooth variation of the transition 
rates across the continuum with distinct endpoint singu- 
larities. The infrared divergence is strong, S^/!{_(n,uj) ~ 
w i/r,-2 ; l j v _ 2 = -1.346 ... as predicted @,§. At the 
upper band edge, the transition rates approach zero in 
what looks like a linear trend. The density D^(n,oj) of 
ipip states is again flat up to near the upper band edge, 
where it has a square-root divergence [Fig. |(c)] . Unlike 
in the case of the tpip* excitations discussed previously, 
here the divergence of the density of states coincides with 
a zero in the transition rates. In the resulting lineshape, 
shown in Fig. ||(d), the divergence is thus suppressed and 
converted into a cusp. 

In the limit of zero external magnetic field (h — * 0), 
the tptp continuum turns into the more familiar two-spinon 

continuum. The exact two-spinon lineshape of S y(n,uj) 

as previously obtained for N = oo via algebraic analysis 
H does indeed exhibit features very similar to those of 
Sli(n,u)) observed here for the first time. In the zero- 
field case, however, the power-law singularities are accom- 
panied by logarithmic corrections. 
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Fig. 3. (a) Continuum of ifiip excitations, (b) Scaled transition 
rates between the ground state and the ijjrp states at q = ir for 
N = 12,16,... ,28 (open circles) and N = 1536 (solid line). 
The dashed line is a fit a + buj 1 ^' 1-2 of the N = 1536 data at 
uj < 0.25. (c) Density of ipip states at q = n for N = 1536. (d) 
Lineshape at q = n of the iptp contribution to S- + (q,u>). All 
results are for M z = N/4. 





Fig. 4. (a) Energy-momentum relation at M z /N = 192/1536 = 
0.125, 384/1536 = 0.25, 576/1536 = 0.375, 1364/3072 ~ 
0.444, 9980/20000 = 0.499 of the ip* branch that is part of 
the lower boundary of the ipip* continuum, (b) Unsealed tran- 
sition rates for the states shown in (a) and those at M z /N = 
9990/20000 = 0.4995, 9995/20000 = 0.49975. 



Our concluding application pertains to the perpendic- 
ular spin fluctuations as described by {q,uj), for which 

we again use the transition rates (£2|). Whereas the func- 
tion S !_((?, uj) dominates the perpendicular spin fluctua- 
tions in weak magnetic fields, it is the function S-\ (q, uj) 

that carries most of the spectral weight in strong fields. 

A set of dynamically dominant collective excitations 
is identified which consists of unbound pairs of ip and ijj* 
quasiparticles as already encountered in the parallel spin 
fluctuations. However, since the perpendicular and par- 
allel spin fluctuation operators reach these excitations in 
different invariant Hilbert subspaces, the spectral bound- 
aries of the iptp* spectrum in (q, uj) are related to those 

shown in Fig. |l](a) by reflection at the line q = tt/2 pq,p|. 

One outstanding feature of the spectral-weight distri- 
bution Stt\q , uj) is the distinct scaling behavior of the 
transition rates for the ip* branch of the lower boundary 
||. Figure ^(a) shows the energy-momentum relation of 
this set of excitations for various values of M z /N. This 
particular branch only exists for M z ^ 0. With M z in- 
creasing from zero, it emerges at uj = 0,q = 0. The 
frequency of its member state at q — increases pro- 
portional to h and the wave number of its member state 
at uj = increases proportional to M z . Upon saturation 
(M z /N — > 1/2), it turns into the branch of one-magnon 
states with dispersion uj(q) = J(l + cosq). In Fi g. (b) 
we plot the unsealed transition rates M^(q) via (|22j ) of 
these states versus q for the same values of M z /N as in 
panel (a). Also shown are two additional values very close 
to saturation. The iV-dependence at fixed M z /N of these 
transition rates is very weak, in particular at small q. We 
have seen that excitations belonging to a continuum have 
transition rates with very different scaling behavior. 

The tp* transition rate at q = is exactly known for 
arbitrary values of M z : 5+_(0, uj) = {2M z /N)2tt6(uj - h). 
It represents a resonant mode of the field-induced magne- 



tization. In the limit h —> hs = 2J, the transition rates 
become independent of q and carry 100% of the intensity, 

S-\ (3, uj) = 2it5(uj — 2 J), which is well understood in the 

context of magnon excitations 0. This trend is very slow 
but clearly visible in the finite- N transition rate data of 
Fig. |(b). _ 

What is perhaps most surprising is that any one of 
these modes carries a nonzero fraction of the spectral 
weight already far below saturation. Hence they manifest 
themselves in sharp resonance lines, separate from the ad- 
jacent continuous spectral weight distributions, with in- 
tensities that become stronger as the magnetic field in- 



creases. 
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